Recovering General Relativity from massive gravity 
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We obtain static, spherically symmetric, and asymptotically flat numerical solutions of mas- 
sive gravity with a source. Those solutions show, for the first time explicitly, a recovery of the 
Schwarzschild solution of General Relativity via the so-called Vainshtein mechanism. 
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It is now well established that the expansion of the Uni- 
verse is accelerating. This can be explained by a non van- 
ishing cosmological constant, raising various questions at 
the fundamental level, or by a ubiquitous, but mysteri- 
ous, dark energy. Instead of introducing such new dark 
components in the Universe, so far only detected via their 
gravitational effects, a logical alternative is to consider 
possible large distance modifications of gravity. His- 
torically, the introduction of the cosmological constant 
has been linked by Einstein himself to the possibility of 
modifying gravity at large distance via a Yukawa decay 
[H, i.e. to give a mass to the graviton. Interestingly, 
"massive gravity" has recently attracted a lot of inter- 
est in this context via constructions relying on models 
with extra-dimensions like the Dvali-Gabadadze-Porrati 
(DGP) model [2J. This model, where gravity is mediated 
via a resonance of massive gravitons, indeed produces a 
late time acceleration of the expansion of the Universe 
via a modification of gravity at cosmological distances 
Q. The consistency of the DGP model has been ques- 
tioned and variants have been proposed 0] ■ On the other 
hand, the way well tested predictions of General Relativ- 
ity (GR) are recovered in models with massive gravitons, 
is highly non trivial. Indeed, a consistent free massive 
graviton is uniquely defined by the Pauli-Fierz theory 
5], which is well know to lead to different physical pre- 
dictions from those of GR, irrespectively of how small 
the graviton mass is. This is the celebrated van Dam- 
Veltman-Zakharov (vDVZ) discontinuity [f|. This, based 
on observations such as light bending around the Sun, 
would be enough to reject any theory of massive gravity 
provided it is well approximated by the Pauli-Fierz the- 
ory. However, soon after it was discovered, it was pointed 
out by A. Vainshtein that the vDVZ discontinuity could 
be cured in a suitable interacting extension of Pauli-Fierz 
theory 0. The purpose of this letter is to show that 
this "Vainshtein mechanism" is indeed working for static 
spherically symmetric solutions appropriate to describe 
the metric outside spherically symmetric sources. Al- 
though a cosmological version of the Vainshtein recovery 
of GR was shown to be present for exact solutions of 
the DGP model [8], nothing analogous was known so far 
for static solutions; e.g. static spherically symmetric so- 
lutions of DGP model are only known locally in some 



approximation scheme Q . Our results are not only of in- 
terest for the specific interacting massive gravity studied 
here, but have also implications for more sophisticated 
models such as the DGP model or its variants [2j, |J] , the 
"degravitation" models [l(| , the Galileon and k-Mouflage 
models [Tljj , as well as possibly the Hofava-Lifshitz the- 
ory Q (see e.g. [3]). 

The theory we will consider has one dynamical met- 
ric g^u non derivatively coupled to a background non 
dynamical flat metric f^ v via a scalar interaction term 
Vmt [/, g] ■ It is defined by the following action 
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where R( g ) is the Ricci scalar built with g^. The in- 
teraction Vi n t[f,g] is chosen such that, when expanded 
at quadratic order around a flat background for g^, one 
recovers the Pauli-Fierz mass term with a graviton mass 
given b y m . There are infinitely many possible such 
choices [14j ; we shall investigate here one of the simplest 
forms, also considered in [3 Eil> given by 

V int = { 9liu - f^)(g aT - f aT ) (g^g vr - gTg") . (2) 

We also consider that matter is minimally coupled to the 
dynamical metric g^ v . The equations of motion, obtained 
varying action ([1]) with respect to the metric <? u „, read 
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where denotes the Einstein tensor computed with 
the metric g M „, is the energy momentum tensor of 
matter fields, and T^fJ is an effective energy momentum 
tensor coming from the variation with respect to the met- 
ric g^v of the interaction term \J—g x Vint ■ 

An ansatz for g^ v appropriate to describe static spher- 
ically symmetric solutions reads 
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where R is parameterizing the distance to the source. 
With such a choice, the solution for g^ v can easily be 
compared to the Schwarzschild solution of GR given by 
Eq. flD with v = v GR = ln(l - Rs/R) and A = 
Xgr = —vqr, where R$ is the Schwarzschild radius of 
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the source. For simplicity, the metric f^ v is taken to 
cover a Minkowski space-time. However one can see that 
when one chooses the ansatz (Q| for g^, one should keep 
one unknown function (i(R) in /„„ in order for the field 
equations not to be overconstrained. Hence, we take f^ v 
in the non canonical flat-space form @, EB, 17 1 



expansion with similar properties was also found for fi). 
For the theory (Q~|[2]) this expansion reads (for R 3> R$) 
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where here and henceforth a prime denotes a derivative 
with respect to R. Note that fj, cannot be measured di- 
rectly because it only appears in and hence, in the 
following we will focus on the physical components {v, A}. 
To describe a static spherically symmetric source, we con- 
sider a matter energy momentum tensor in the perfect 
fluid form reading T" — dia,g(—p,P,P,P), where the 
pressure P and energy density p depend only on R. We 
also assume that this energy momentum tensor vanishes 
outside the radius Rq of the source. For such a perfect 
fluid and the above ansatz Q), the conservation equation 
reads 



P' = -j(P + P ). 



(6) 



We look for asymptotically flat and everywhere non sin- 
gular solutions. In particular, we impose boundary con- 
ditions such that /j,, and v vanish at infinity as well as A 
and // in R = 0. 

If one linearizes the vacuum equations of motion, one 
finds that v and A are respectively given by 
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where C is some integration constant. Notice that for 
R <C mT 1 one gets that vl ~ — 2Al irrespectively of the 
graviton mass. This is in contrast with the corresponding 
relation for GR stating that at distances larger than Rq 
one has vqr. = —\gr, and is just another way to phrase 
the vDVZ discontinuity. The constant C should be fixed 
by including a source. 

Following Vainshtein, one can then show that the solu- 
tion {vli^l} obtained by linearization ceases to be valid 
at a distance smaller than the Vainshtein radius, Ry, 

given by Rv — (m~ 4 i?s) For a graviton with cosmo- 
logically large Compton length m _1 , the Vainshtein ra- 
dius of the Sun is much larger than the solar system size 
impeding the use of the linearized solution to compare 
with solar system observables. Vainshtein argued that, 
below Ry, a recovery of the standard solution of GR 
was possible. Indeed, he found a local expansion of the 
functions {A, v\ around their GR expressions {Xgr, v gr\ 
that stays valid at small m and for distances R <C Ry (an 



v ~ v GR 
A — \gr 



(rni?)V(8i?s)/(81i?) + 0(m 4 ), (9) 
-(mR) V(2i?s)/(9i?) + 0{m A ). (10) 



In contrast, the linearized expressions {Xl^l} are 
nowhere valid when m goes to zero, because Rv is in- 
finite in this limit. However it was pointed out by Boul- 
ware and Deser [3| that there was no warranty that the 
two expansions found by Vainshtein (for R ^> Ry and 
R <C Rv) could be joined into a non singular solution of 
the field equations. This question remained unanswered 
for a while and it was more recently investigated in Ref. 
[l6| which concluded, on the basis of a numerical investi- 
gation of the field equations, that no such solution existed 
and hence that the Vainshtein mechanism, i.e. the non 
perturbative recovery of solution of massless gravity, does 
fail in the massive gravities defined by the class of action 
(JT|) , and in particular for theory ([TJ2]). Here we reexam- 
ine this question in the light of a recent work (l9l | where 
we considered static spherically symmetric solutions of 
the field equations ([3]) in the the so-called decoupling 
limit (DL). Note that the theories in the class (UJ) suffer 
from the Boulware-Deser instability [l|| associated with 
a mode that starts propagating at the nonlinear level. 
This is probably enough to make those theories unreal- 
istic, but is not really worrisome for the present work, 
since we simply view the theory ([T][2]) as a tool to inves- 
tigate the validity of the Vainshtein mechanism, and do 
not advocate for a realistic application of this particular 
theory. 

The DL enables to focus on the nonlincaritics respon- 
sible for the crossover behaviour arising at the Vain- 
shtein radius Ry and is obtained by sending together 
Mp and m _1 to oo while keeping the energy scale A = 
(m 4 Mp) 1 / 5 , as well as Rs x Mp, fixed. The scale A is as- 
sociated with the strongest self-interaction in the scalar 
sector of the model [l5| . The DL should provide a good 
description of the solutions (provided those exist) in the 
range of distances Rs <C R -C m -1 . Notice that this 
range of distance covers in particular the crossover dis- 
tance Ry which stays constant in the DL. In [3] we 
obtained numerical solutions of the DL field equations 
which are showing a transition between a R <C Ry regime 
recovering linearized GR behaviour and a R Ry 
regime, where the leading behaviour of the solutions is 
given by Eqs. ([7][8]) where one sends my. Rio zero. Such 
solutions were shown to be everywhere non singular and 
to exist only for a class of theories of the kind |T]) includ- 
ing the one defined by Eq. ([2]). Note that, in the DL 
the nonlinearities important in GR at the Schwarzschild 
radius are absent as well as the terms responsible for the 
Yukawa decay at large distances. A crucial property of 
the DL, is that, in vacuum, we found infinitely many 
solutions decreasing at infinity. This was shown both 
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FIG. 1: Plot of the metric functions — v and A vs. R/Rv, 
in the full nonlinear system and the decoupling limit (DL), 
with a star of radius _Rq = 10~ 2 Ry and m x Rv = 10~ 2 . For 
R <C Rv, the numerical solution is close to the GR solution 
(where in particular v ~ — A ~ —Rs/R for R > -Rq). For 
R 3S> Rv, the solution enters a linear regime. Between Rv 
and m _1 , where the DL is still a good approximation, one 
has v ~ — 2A ~ —4/3 x Rs/R. At distances larger than mT 1 
the metric functions decay a la Yukawa as appearing more 
clearly in the insert. The latter shows the same solution but 
for larger values of R/Rv, and in the range of distance plotted 
there, the numerical solutions are indistinguishable from the 
analytic solutions of the linearized field equations Eqs. (|7I8[) . 



by numerical integration and by using Ecalle's theory of 
resurgent functions and, in particular, allows one to ac- 
commodate the presence of a non vanishing source. 

The existence of non singular solutions in the DL is 
not enough to conclude on the validity of the Vainshtein 
mechanism in the full theory. First because those solu- 
tions are only covering at best a limited range of distances 
Rs <C R <C mT 1 (this interval being sent to [0, +oo] tak- 
ing the DL) and second because it could be that non- 
linearities present in the original field equations, but left 
aside by the DL, destabilize the DL solutions enough to 
make them become singular. 

In this letter we report results of the numerical inte- 
gration of the full nonlinear field equations (J3j) — dOJ) for 
the theory fl][2|) and metric ansatze (@][5|). Those equa- 
tions, whose explicit (but rather long) form can be found 
in Ref . [l6| . reduce to a set of quasilinear ordinary dif- 
ferential equations of first order in A and v and of sec- 
ond order in /i. We included a static spherical source 
of radius Rq and constant density pq. Hence, defin- 
ing the Schwarzschild radius Rs of the source as Rs = 
Mp 2 RqPq/3, there are only two dimensionless parame- 
ters in the model which can be taken to be Rq / Rv and 
a = m x Ry . We were able to find everywhere non singu- 
lar asymptotically flat solutions, {v, A, fi, P}, by solving 
numerically the field equations, using both shooting and 
relaxation methods. To our knowledge, the existence of 
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FIG. 2: Plot of the ratio of v to —A vs. R/Rv, with a star 
of radius Rq/Rv = 10~ 3 and m x Rv = 10 -3 . This shows 
the transition, at the Vainshtein radius, between a GR regime 
where v ~ —A to a regime where v ~ — 2A. At larger distance, 
the solution features the expected Yukawa cutoff. 



such solutions was unknown and shows for the first time 
explicitly that the Vainshtein mechanism can work in a 
theory of massive gravity. The shooting approach, which 
is based on a direct Runge-Kutta integration, is quite 
difficult to implement since the initial conditions need to 
be extremely finely tuned in order to obtain non singular 
solutions. This was already the case in the DL and 
could explain the findings of Ref. [l^|. We have been 
able to integrate the system of equations outward from 
a point deep inside the source up to a point further than 
the Vainshtein radius (typically R ~ 3i?y). We have also 
been able to integrate inward starting from a point far 
from the source (e.g. R ~ 3Rv) up to the core of the 
source. In fact the integration of the vacuum equations 
shows that, as in the DL, one is able to obtain infinitely 
many decaying solution by slight changes of the bound- 
ary condition. In all the integrations, we have made an 
extensive use of local analytic series expansions of the 
solution in order to fix the initial conditions as precisely 
as possible. Even in cases where we could check that 
those expansions were only asymptotic expansions (i.e. 
that the series were not convergent), they were very use- 
ful to set with enough precision the boundary data of the 
numerical integrations. Details of the numerical integra- 
tion will be reported elsewhere [2(| ■ In order to be able 
to obtain a wide enough range of integration (typically, R 
runs between R = 10~ 5 i?y and R ~ 3Ry) the required 
relative precision on the initial conditions is usually of 
order 10 -8 and often even more stringent. These results, 
obtained through a Runge-Kutta integration, have been 
confirmed by the use of a relaxation method. In particu- 
lar the relaxation approach allowed us to solve the system 
on a wide range of parameters and distances (typically 
for R between and ~ 100 Rv). Both methods perfectly 
agree, demonstrating the robustness of our numerical in- 
vestigation. Our main finding is that for a € [10 _3 ,0.6] 
and sources that are not too compact (i.e. Rq > 5Rs), 
a non singular solution can be found, which is asymptot- 
ically flat and has the right boundary conditions at the 
origin R = 0. It features a recovery of GR at distances 
R <C Rv, where we also checked numerically that the 
first correction to the GR behaviour was given by the 
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FIG. 3: Plot of — vx a~ 4 vs. R/Rv (the a -4 factor is included 
for convenience such that, in the decoupling limit (DL), all 
plotted theories would exactly coincide) for three different 
values of a = m x Rv and source of radius Rq = 10~ 3 Rv- 
The solution with a = 0.005 corresponds to a spherical source 
of size and density close to that of the Milky-Way and a gravi- 
ton Compton length of the order of the Hubble radius. In the 
range of distances plotted, this solution (with a = 0.005) is 
well approximated by the DL solution. The insert shows a 
zoom on small distances. There, the behaviour of our solu- 
tions agrees with the one of GR, and for a — 0.1, a value which 
belongs to a parameter range investigated in Ref. [l6| ]. one 
can see that the solution departs from the DL, emphasizing 
the role of the non-linearities of GR which become important. 



form given in Eqs. ([gi fTU)) . It is also well approximated 
by the solution obtained in the DL in the expected range 
of distances Rs <C R -€1 m _1 . In particular, at distances 
above Ry (but below m _1 ) the behaviour of the solution 
is given by the regime 



-2A - -4/3 x R s /R, 



(11) 



as obtained from Eqs. ([7][5]). Note that the constant 
C is found to be ~ 4/3 x Rs by the numerical integra- 
tion. This value of C agrees with that obtained in the 
linearized theory for a non relativistic point-like source 
placed at the origin, and can be explained by the presence 
of a scalar polarization also responsible for the vDVZ dis- 
continuity. At distances above m _1 the solution features 
the exponential Yukawa falloff. We checked numerically 
that this falloff agrees with the one given in Eqs. (jTJIHI) 
as well as the fact that our solutions were continuously 
connected to flat space-time, when the mass of the source 
is decreased either by keeping the density fixed but de- 
creasing the radius or by keeping the radius fixed but de- 
creasing the density. Those results are depicted in Figs. 

To summarize, we found static, spherically symmetric, 
and asymptotically flat numerical solutions of massive 
gravity showing, for the first time explicitly, a recovery of 



the usual Schwarzschild solution of General Relativity via 
the so-called Vainshtein mechanism. Some issues deserve 
more investigations. These include the classical stability 
of our solutions as well as an understanding of highly 
relativistic objects and Black Holes. 
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